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Aims. To clieck whetlier the polar angle distribution of QSOs around nearby spiral galaxies is isotropic or not. 
^ Methods. A statistical analysis of the polar angle distribution of large samples of QSOs from the SDSS survey and Monte Carlo simulations to 
calculate their significance are carried out. 
j Results. There is a clear excess of QSOs near the minor axis with respect to the major axis of nearby edge-on spiral galaxies, significant at a 
Q,^ level 3.5cr up to angular distances of ~ 3° (or ~ 1.7 Mpc) from the centre of each galaxy. The significance is increased to 3.9cr with the z > 0.5 
^— s ■ QSOs, and it reaches 4.8cr if we include galaxies whose circles of radius 3 degrees are covered by the SDSS in more than 98% (instead of 
\^ [ 100%) of the area. 

Conclusions. Gravitational lensing in the halo of nearby galaxies or extinction seem insufficient to explain the observed anisotropic distribution 
of QSOs. The anisotropic distribution agrees qualitatively with the predictions of Arp's models, which claim that QSOs are ejected by galaxies 
along the rotation axis, although Arp's prediction give a distance of the QSOs ~ 3 times smaller than that found here. In any case, a chance 
fluctuation, although highly improbable, might be a possibility rather than a true anisotropy, and the present results should be corroborated by 
O ' other groups and samples, so we prefer to consider it as just a first tentative detection. 

^ Key words, quasars: general - Galaxies: statistics - Catalogs - distance scale - Gravitational lensing 
■ 

> : 

k>( 1 . Introduction The recent analysis carried out by Scranton et al. (2005), who 

^ used Sloan Digital Sky Survey (SDSS) photometric data con- 

d The first hints of a possible relationship between nearby galax- j^j^j^g 2 x 10^ quasars and ~ 1.3 x 10« galaxies, showed 

les and high redshift QSOs came from Arp (1966, 1967), who ^^^^ ^j^^ amplitude and sign of the angular coiTelation func- 

observed radio sources across active galaxies that were iden- ^j^^ heWeen both sets of objects depend on the magnitude 

tified as QSOs. Since the pioneering work by Burbidge et al. ^-^^ considered for the sample of QSOs. From this correla- 

(1971), several groups have demonstrated the existence of an tion, Scranton et al. (2005) proposed an ad hoc halo distribution 

angular correlation between samples of QSOs and low redshift function compatible with a cross-correlation of very small am- 

galaxies (e.g. the reviews of Guimaraes 2005 and Burbidge pjij^^g (^^^ < 0.04) of faint galaxies with QSO candidates 

2001 respectively of orthodox and heterodox approaches). This selected photometrically (5% of this sample are not QSOs; 

correlation extends to angular scales of ~ 1 degree. Although Richards et al. 2004). It is small because the mean separa- 

weak gravitational lensing by dark matter has been proposed ^-^^ ^^^^^ g^j^^i^^ ^^^^jj ^^^-^-^^ con-elation of 

to be the cause of these coiTelations, many authors have found q^q^ ^^^^^^ ^ g^j^^y ^-^^^^^ ^-^^^ ^j^^ contamination of 

this hypothesis insufficient to explain the correlations (Benitez ^^^^ ^^^^^^ qSq^ belonging to other galaxies. This still seems 

et al. 2001 ; Gaztanaga 2003; Nollenberg & Williams 2005). insufficient to solve the most important problem of the correla- 

tion found between QSOs and the nearest and brightest galax- 
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ies (Kovner 1989); for instance, for the high amplitude angular 
correlation found by Chu et al. (1984) wgq ~ 5. Although in- 
completeness could be responsible for these correlations, it is a 
matter that requires further study. 

The existence of a correlation between samples of objects 
with different redshifts has been advocated by the supporters of 
non-cosmological redshifts (e.g. Burbidge 1999) as clear evi- 
dence of physical association and as proof that QSOs are being 
ejected by low redshift galaxies. One variant of the model (e.g. 
Arp 1998a, ch. 3; Arp 1999a) assumes that QSOs are ejected 
along the rotation axis of the parent galaxy, decreasing in red- 
shift as they move outwards and eventually becoming normal 
galaxies. 

There have been claims of anisotropy in the distribution 
of QSOs around nearby galaxies. For instance, there are con- 
figurations of QSOs aligned along the minor axis of a central 
Seyfert whose probability of being accidental is 10"^ to 10"'° 
(Arp 1998b, 1999a). Arp & Hazai'd (1980) have reported three 
QSOs in a straight line, together with another set of three QSOs 
in a straight line, a configuration that is very difficult to explain 
as random. Arp (1999b) finds that the probability of having 
six out of six QSOs aligned within ±15° of the minor axis of 
NGC 5985 to be only 10 ^* to 10 ''. Arp & Russell (2001) found 
that the bright radio quasars 3C37 and 3C39 are paired across 
the centroid of the disturbed galaxy pair NGC 470/474 in an 
arrangement with a probability of only 2 x 10"^ of being acci- 
dental. There is anisotropy in the radio QSO distribution at high 
flux densities (Shastri & Gopal 1983): the number of QSOs on 
one side of the M33 region is far greater (~ 1 Icr) than that of 
the diametrically opposite region. This alignment of QSOs with 
the parent galaxy would involve conservation of momentum in 
their ejection in the heterodox Arp/Burbidge/Narlikar scenario 
(Narlikar & Das 1980). 

In this paper we undertake the first study of the polar angle 
distribution of QSOs around bright nearby galaxies in statis- 
tical terms using large samples of both types of objects. With 
the release of the SDSS survey, a systematic search for this ef- 
fect is now possible. The motivation of this research is to test 
one of the predictions given by the above authors concerning 
the ejection model. The predictions are as follows: if the QSOs 
were background objects, no preferential orientation in princi- 
ple should be found; however, if the QSOs were ejected by the 
galaxies, as Arp proposed (or alternatively a gravitational ef- 
fect were predominant in the rotation axis of a galaxy), a higher 
concentration of objects near the rotational axis of the galaxy 
would be expected. 

Section 2 presents the data and sample selection; the anal- 
ysis and statistical tests are presented in Sections 3 and 4, 



Section 5 discusses some details of the method to measure the 
anisotropy, and in Section 6 we interpret the results. 

2. Data and sample selection 

Two sources of public data have been used in this work: 

- Third Reference Catalog of Bright Galaxies (de 
Vaucouleurs et al. 1991, updated on 16 February 1995; 
hereafter RC3): this catalogue is complete for galaxies 
with apparent diameters greater than 1 arcmin at the D2s 
isophotal level and total Z?-band magnitudes Bj < 15.5, 
with a redshift z < 0.05. It contains a list of 23 01 1 galaxies 
over the whole sky. 

- Spectroscopic catalogue of the Third data release (DR3) 
of the Sloan Digital Sky Survey (Abazajian et al. 2005; 
hereafter SDSS): this catalogue covers 4188 square de- 
grees and contains 51 027 QSOs. The spectra have signal- 
to-noise ratios SNR > 4 per pixel at nig = 20.2. From 
these QSOs, we took the vetted subsample by Schneider 
et al. (2005) in which only those objects with M,- < -22 
are selected (this implies that the low redshift QSOs are 
removed, and that the redshift range is from 0.08 to 5.41; 
their redshift is therefore always higher than those of the 
RC3 galaxies) and many of them have been checked manu- 
ally to be QSOs (or otherwise removed). The total number 
of QSOs of this vetted catalogue is 46420. We prefer to 
use the spectroscopic data instead of the photometric data 
used, for example, by Scranton et al. (2005) with a colour 
selection technique to separate the QSOs for reliable iden- 
tification and an accurate estimate of the redshift. While we 
were working with these data, a fourth release of the survey 
was delivered that covered an area 14% larger (Adelman- 
McCarthy et al. 2006), but an updated vetted subsample has 
not been produced yet (and will not be produced until a fifth 
release of the data; D. R Schneider, priv. comm.). In ^we 
show some results of the unvetted 4th release. 

We are aware that the DR3 Quasar spectroscopic catalogue 
is not a statistical sample. Its QSOs do not all follow the same 
target selection, and the completeness depends on the positions, 
magnitudes, redshifts, etc., of the QSOs. Nonetheless, it is ap- 
propriate for the statistical analysis presented here because that 
completeness is related to the intrinsic properties of the QSOs 
and not to the orientation of nearby galaxies. Possible bias due 
to incompleteness and other systematic effects in the sample 
will be considered and quantified through Monte Carlo simula- 
tions in Section^] 

We take the RC3 catalogue and select all spiral galaxies 
with information on radial velocity (v3k), the angular radius of 
the 25th mag isophote (^25), the ratio between the major and 
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minor axes and position angle, and whose coordinates are such 
that a circle of radius equal to ^max around each galaxy was to- 
tally covered by the 3rd release of the SDSS survey (see ^5.31 
for the use of circles partially covered). Spiral galaxies were 
chosen because the inclination of their discs is easily derived 
from the projected axial ratio (cos / = ^). We reject the very 
nearby galaxies with V3k < 700 km/s (because their distances 
cannot be determined from their redshifts and because they 
would also extend too much in angular size with respect to the 
average sample; in any case, they are very few in number, and 
their inclusion would not significantly alter the statistics in this 
paper). We also select edge-on galaxies within the range ±25° 
(sin / > 0.906). There is no convenient a posteriori selection of 
galaxies: the selection of edge-on galaxies is necessary because 
they allow a clear separation between the rotation axis and the 
disc of the galaxy (the projection of the rotation axis is perpen- 
dicular to the disc of the galaxy and has a negligible component 
in the line of sight), and only spiral galaxies provide informa- 
tion for deriving the inclination to make possible the separa- 
tion of edge-on galaxies. A range of 25 degrees (/ > 65 deg) 
was taken in order to achieve a high enough number of galaxies 
within the nearly edge-on criterion (sin/ > 0.906). If the lim- 
iting inclination were slightly larger or smaller, the statistical 
results given in this paper would be similar. If we selected only 
galaxies with an inclination very close to 90 degrees, we would 
have very few galaxies in our sample and the statistics would 
be much poorer; if we took the limiting inclination much lower 
than 65 degrees, we would introduce noise from nearly face-on 
galaxies. 

We select the 44975 QSOs with nig < 21 within an an- 
gle of 025 < < flmax from cach galaxy (the lower limit in 
distance avoids possible effects of galactic extinction in back- 
ground QSOs). There is no overlap in the redshift of the galax- 
ies (z < 0.05) and QSOs (z >~ 0.08). Figure [l] shows the posi- 
tion of the resulting samples for the case ^^ax = 3°. 



3. Analysis 

We compute the position angle of QSOs with respect to the 
minor axis of each galaxy. For ^^ax = 3° the total number 
of QSO-galaxy pairs is 25 176 for 71 RC3-galaxies. These 71 
galaxies follow these constraints, with an average distance of 
46 Mpc (derived from V3k, assuming a Hubble constant of 72 
km s ' Mpc ') and a median distance of 32 Mpc. Some QSOs 
are counted more than once because they are within the radius 
^max of two or morc galaxies (see Fig.^, but their relative po- 
sition angles with respect to the different galaxies with which 
they are associated are uncorrected, so they count as indepen- 
dent measures (see discussion in ^5.2l i. 
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Fig. 1. Position on the sky of the sample of objects analysed 
in this article. Large open circles represent galaxies, and small 
dots QSOs. The plot corresponds to the case ^^^ax - 3°. 
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Fig. 2. Histogram representing the number of QSOs with nig < 
2 1 within a circle of radius equal to 3 degrees of 7 1 RC3 galax- 
ies as a function of position angle with respect to the minor 
axes of these galaxies. 
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Fig. 3. a vs. Omax for QSOs with nig < 21 and any redshift (the 
points are not independent). 

A histogram of all counts vs. relative position angle with 
respect to the minor axis of the galaxy is shown in Fig. |2 for 
^max = 3°. Clearly, there is a decrease in counts as we move 
from the minor to the major axis. The excess of counts in the 
minor axis direction is a; 13% higher than in the major axis 
direction. A fit for the counts per bin vs. position angle (PA) 
with a linear function 



C = C()(l +ar), 



(1) 



/ PA \ 

where r = |^1 - j (between and 1) 



gives Co = 1311 ± 13, a = 0.132 + 0.017. This means a de- 
tection of anisotropy at the 7.8cr level, although this is only the 
statistical significance of the survey, which does not represent 
the statistical significance of the real (unbiased) distribution of 
QSOs, as discussed in Section^ 

Other values of fl^ax (see Fig. |3} show the existence of 
anisotropy but the best significance is for O^a^ ^ 3°. For lower 
angles, there seems to be a lower value of a, although for 
^max < 1° the number of QSOs is too low to draw any con- 
clusions. For 0max > 4° the number of galaxies with total SDSS 
covered area in the circle with radius flmax is much lower (for 
^max = 4.5° we have only 12 galaxies). In any case, the detec- 
tion ofai^O (anisotropy) is clear, and the projected linear scale 
of around 1-1.5 Mpc (around 1.5-2.0 Mpc without projection) 
seems to have the highest ratio of anisotropy. 

The excess in the range 0^5° (13 025 cases) over the range 
45-90° (12 151 cases) in Fig.Qis £45 = -h874 cases, a ratio of 
excess of 7%. These extra 874 sources are preferably placed in 



the magnitude range nig > 19.4 and the range 0.2 < logwil + 
z) < 0.5 (i.e. 0.6 < z < 2.2), as can be observed in Figures |4] 
and Is] Curiously, the QSOs with z <~ 0.5 significantly show 
the opposite trend on average: the excess is towards the major 
axis. The peaks with significant excess (higher signal/noise) or 



relative excess are more or less for lUg 
respectively. 



20.1 and z ~ 0.6 



Figure |6l also shows how the QSOs with anisotropy are 
preferentially those with nig > 19.4. The limit is near the max- 
imum of the differential QSO counts; over 19.2 the counts are 
lower because of appreciable incompleteness (otherwise the 
differential counts should increase monotonically). One might 
suspect that the anisotropy has something to do with the in- 
completeness/bias of QSOs; however, it will be demonstrated 
in ^4.11 that the main cause of the anisotropy is not that; at 
least the effects of incompleteness/bias, if there are any, are 
not enough to explain the observed anisotropy. Examining Fig. 
|6l one should also realize that there are two structures in the 
counts: two overlapping peaks, one with a maximum at nig » 
19.2 and another at nig » 20.2. Apparently, it is the second 
group of QSOs that is responsible for the anisotropy, and this 
is shown over irig > 19 A because this is the range where the 
number of QSOs in the second group is relatively significant. 
We suggest this to be the main reason for the manifestation 
of the anisotropy over nig > 19.4 rather than the incomplete- 
ness. This second peak corresponds mainly to the criterion of 
target selection in the SDSS survey (Schneider et al. 2005) of 
serendipity sources with unusual colours; other targets do not 
make a such an important contribution (see Fig.0. Values of a 
with the different target selections are given in Table ^ When 
we take the QSOs that appear in the Serendipity target and ex- 
clude the other targets, we get a value of ff = 0.380 ± 0.055, an 
excess of 38% towards the minor axis instead of the average of 
13% with the whole sample. 

If we do a new analysis of the anisotropy with the added 
constraint that nig > 19.4, 0.2 < logu)(l + z) < 0.5 (a total 
of 13 308 cases), we get 7509 cases of QSO-galaxy associa- 
tion, and a to 0.313+0.042. This represents more than a 30% 
excess towards the minor axis with respect to the major axis. 
The analysis of the anisotropy only with the further added con- 
straint that nig > 19.4 gives 10204 cases of QSO-galaxy asso- 
ciation with a - 0.237 ± 0.036. If we restrict the analysis with 
the added constraint that z > 0.5 (log(l + z) > 0.176) we get 
22 953 cases of QSO-galaxy association and a = 0.156+0.017. 
Although this simple constraint (in which we remove only 9% 
of the associations with respect to the analysis of the full sam- 
ple) does not give a much higher value of a; it is much more 
significant statistically, as we comment in Section|4] The results 
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Table 1. Anisotropy of QSOs with different constraints 



^max 


Other constraints 


Number of pairs 


a 




0.5° 




2807 


-0.011 H 


t 0.056 


1.0° 




8698 


0.004 ± 


0.049 


1.5° 




14902 


0.036 ± 


0.039 


2.0° 




20562 


0.067 ± 


0.026 


2.5° 




23819 


0.123 ± 


0.017 


3° 




25176 


0.132 ± 


0.017 


3° 


nig > 19.4 


10204 


0.237 ± 


0.036 


3° 


ra^ > 19.4, 0.2 < log,o(l + z) < 0.5 


7509 


0.313 ± 


0.042 


3° 


z>0.5 


22953 


0.156 ± 


0.017 


3° 


Low-z QSOs target 


15434 


0.064 ± 


0.026 


3° 


High-z QSOs target 


5055 


0.087 ± 


0.044 


3° 


FIRST (radio) target 


1184 


0.203 ± 


0.085 


3° 


ROSAT (X-ray) target 


1728 


0.031 ± 


0.087 


3° 


Serendipity target 


14445 


0.179 ± 


0.030 


3° 


Serendipity (only) target 


5859 


0.380 ± 


0.055 


3° 


Star target 


304 


0.164 ± 


0.200 


3° 


Galaxy target 


205 


-0.303 H 


t 0.207 


3.5° 




21286 


0.102 ± 


0.033 


4° 




14130 


0.059 ± 


0.032 


4.5° 




9328 


0.148 ± 


0.047 




Fig. 4. Ratio of relative excess of QSOs in the range of position 
angles with respect to the minor axis 0^5° over the range 45- 
90° versus magnitude of QSOs. 




Fig. 5. Ratio of relative excess of QSOs in the range of position 
angles with respect to the minor axis 0^5° over the range 45- 
90° versus redshift of the QSOs. 



of these and other analyses with different ^max are summarized 
in TableC] 

The two-dimensional plot of the polar and radial distribu- 
tion for two of the cases analysed are presented in Figure |8] 
Apparently, QSOs are not randomly distributed around galax- 



ies but show a quite significant trend to follow the minor axes 
of those galaxies. Since the median distance of the galaxies is 
32 Mpc, the median linear projected maximum distance corre- 
sponding to flmax - 3° is 1.7 Mpc. 
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Fig. 6. Differential QSO counts in bins of 0.2 mag around the 
71 galaxies that follow flmax = 3°. The angle is with respect to 
the minor axis. 
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Fig. 7. Differential QSO counts in bins of 0.2 mag around the 
71 galaxies which follow ^^ax = 3° with an angle respect to 
the minor axis of less than 45 degrees. Different SDSS samples 
attending to different criteria to preselect the targets are used. 



Figure|9lpresents for the 71 galaxies considered in the case 
^max - 3° the number of QSOs at ±45° from the minor and 
major axes respectively. The values are also given in Table |2] 
There is in Fig.|9]a different degree of dispersion of anisotropics 
depending on the number of galaxies (galaxies with number 
less than 20 and higher than 60 present in general more disper- 




Fig.8. Top : Counts of QSOs (total: 25 176) as a function of 
position with respect to the minor/major axis angles of the cor- 
responding 71 galaxies (see text). Counts were plotted in bins 
of 0.2° X 0.2° (with an average of around 140 QSOs per bin) 
and smoothed/interpolated. The clearer colours indicate higher 
density. Note that towards the minor axis the average density 
of QSOs is greater than towards the major axis. Most of the 
fluctuations are presumably statistical (Poissonian). Bottom : 
The same but with the extra constraint z > 0.5 (total: 22953 
counts). 

sion in anisotropy). This may be due to the different levels of 
completeness for the different galaxies; the galaxies are num- 
bered in Fig.|9]in right ascension order, so this most probably 
means that there are some regions of the sky in which the com- 
pleteness is lower, and the number of background sources is 
also lower, as observed. 
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Table 2. Distribution of SDSS QSOs with < 21 within 025 < 
< 3° around the 71 selected spiral galaxies (with / > 65° 
that have complete coverage of circles of radius 0„ax = 3°). 
£45: excess of QSOs in the region which are at position angles 
less than 45° from the minor axis over those that have posi- 
tion angles greater than 45° from it (with Poissonian error). 
The galaxy type or environment is taken from the SIMBAD 
database classification. Distance (d) derived from the redshift 
with a Hubble constant 72 km/s/Mpc. PA is the position angle 
of the major axis of the galaxy. 



# PGC-# Type/env. 





POP-* 


y r 1 * • 




PA (°^ 




C45 


1 

1 








79 


497 




9 

z. 


Z.JZ. JZ, 




1 1 1 








J 


zj / 1 




J J .0 


1 90 


j\jy 




4 


25910 


T ow S Rr 


32.6 


98 


329 


-15(+18) 


5 


26232 


Interact. 


27.0 


160 


462 


+58(±21) 


6 


26238 


in Pair 


27.8 


125 


462 


+44(±21) 


7 


26246 


in Pair 


34.4 


48 


323 


-25(+18) 


8 


26351 




9.9 


140 


307 


+25(+18) 


9 


26563 




40.1 


4 


413 


+37(+20) 


10 


26631 


Interact. 


40.5 


50 


299 


+49(±17) 


11 


26856 




46.6 


123 


548 


-58(±23) 


12 


27248 




32.3 


47 


436 


+50(±21) 


13 


27734 


in Pair 


32.8 


72 


460 


-16(±21) 


14 


28010 




30.6 


127 


463 


-21 (±22) 


15 


28033 


in Group 


88.3 


68 


503 


-41 (±22) 


16 


28148 




30.7 


15 


505 


+49(±22) 


17 


28741 




34.5 


18 


482 


+64(±22) 


18 


28939 




33.1 


151 


529 


+57(±23) 


19 


30364 




99.7 


62 


384 


-36(±20) 


20 


30885 




54.3 


111 


299 


+ 17(±17) 


21 


31037 




21.2 


168 


318 


-2(±18) 


22 


31269 




101.3 


156 


323 


-1(±18) 


23 


31608 




96.6 


125 


326 


+ 14(±18) 


24 


32153 




18.6 


159 


313 


+29(±18) 


25 


32564 




42.2 


78 


328 


+32(±18) 



Two of the galaxies show a level of anisotropy higher than 
3<T. The galaxy with highest anisotropy is PGC 36192 [NGC 
3795]. See the distribution of QSOs around it in Figure^] The 
excess £45 (excess of QSOs of this galaxy in the region that 
are at position angles less than 45° from the minor axis over 
those that have position angles larger than 45° from it) is -1-99 
(5.2cr over zero), and for this particular case a - 2.02 + 0.28, 
an excess of 200% towards the minor axis with respect to the 
major axis. 

Among the 71 galaxies for flmax - 3°, there are (accord- 
ing to the SIMBAD database, see Tabled}: one LINER, two 
interacting galaxies, four low surface brightness galaxies, 15 



26 
27 
28 
29 
30 
31 
32 
33 
34 
35 
36 
37 
38 
39 
40 
41 
42 
43 
44 
45 
46 
47 
48 
49 
50 



32570 
32714 
33234 
33766 
33964 
34030 
34971 
35037 
35249 
35675 
35900 
36102 
36192 
38117 
38120 
38188 
39832 
41618 
42255 
42336 
42689 
42791 
42942 
42975 
42998 



in Group 



in Cluster 
in Cluster 
in Cluster 
Low S.Br, 
in Pair 



in Group 
in Group 

LINER 



Cont. Table|2J 

d (Mnc) 


r/\^^ ) 


QSOs 




1 Q 

\.y.j 


1 ^0 


jjj 


-i-9S('-t- 1 


z,o.o 






4- U+ 1 R"! 


70 7 


90 


'^07 


4_^f+ 1 R"! 


?Q 9 




9Q9 


j_ 1 (.(+ 1 71 

^ 1 O^It V I } 


49 It 

Hz.. J 


JO 


9Q7 

Ay 1 


+3\±l 1 ) 


1 9 n 


oU 






ij.y 


1 49 


1 9 








^04 




9S 


1 98 


^no 


-i-Sf'-t- 1 71 
±o^± i- 1 ) 


IQ S 


UU 


4SQ 


1 ^C+91 1 


1 Q 9 


Q 

y 


1 




OU-H- 


1 "^9 


"^OS 




L / .J 


J J 


D 1 J 


4.qQC+ 1 qi 


91 Q 


1 dQ 

L'-ry 


980 


w+ 1 1\ 
vy± 1 1 ) 


87 S 
/ .0 


1 dCi 


9AS 
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galaxies in pairs/groups/clusters, and 49 normal galaxies. We 
looked for signs of correlation of the galaxy Hubble type with 
the anisotropy and found nothing. There is no correlation of 
£45 with the position angle of the galaxy (see Table|5}- There is 
some trend of higher anisotropy with lower distance (see Table 
E): only 5(±5)% (1/19) of the galaxies with > 50 Mpc have 
£45 > 2cr(£'45), while for 27(±7)% (14/52) of galaxies <i < 50 
Mpc with £45 > 2cr(£'4g). This may be because 3 deg is too 
great a distance for galaxies with d > 50 Mpc and the possible 
effect of anisotropy is diluted in them. There also seems to be 
a tendency for interacting galaxies to have high anisotropy (see 
Table|3, but good statistics cannot be carried out with only two 
galaxies out of 7 1 . And we have only one Seyfert (LINER), so 
neither can we comment on statistics of this type. 

4. Statistical significance of the anisotropy 

The spectroscopic survey SDSS is not complete for all redshifts 
up to the magnitude nig < 21. Moreover, we have observed that 
the completeness depends on position in the sky because the 
target selection varies from region to region. There are, for in- 
stance, important gradients of QSO density with the right as- 
cension. Moreover, there are selection effects depending on the 
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redshifts of the sources; SDSS is not equally capable of de- 
tecting sources at all redshifts, etc. Nonetheless, it is appropri- 
ate for our statistical analysis because that completeness has 
nothing to do with the orientation of nearby galaxies; that is, 
whether there are more or fewer QSOs has nothing to do with 
the presence of a galaxy some degrees away and even less with 
the orientation of the latter 

There will be a global effect in the dispersion from the 
isotropy due to the different gradients and systematics, which 
will be higher with respect to a perfect statistical survey, so 
we cannot trust the statistical error bars assuming a random 
distribution. For instance, the 7.8cr level which we derived to 
be the significance of non-zero slope of the fit of expression 
Q applied to Fig. |2l would indicate the level of anisotropy of 
the Schneider et al. (2005) sample, which is both due to the 
anisotropy in the real distribution and to the systematic devia- 
tions of a homogeneous sample. We can also calculate this level 
of significance through Kolmogorov-Smirnov test applied to 
the same data, which also gives a low probability of being com- 
patible with an angle-independent distribution: Pk-s - 2x10"*^ 
(equivalent to 5.6 cr in a Gaussian distribution). Or simply 
we can count in the same data the excess in the range 0^5° 
(13025 cases) over the range 45-90° (12151 cases) in Fig.|2l 
£45 - -1-874 cases (Poissonian error ±159), i.e. 5.5cr. In order 
to distinguish which part of the significance of the anisotropy 



Fig. 9. Top : Number of QSOs at +45° from the minor axis 
(open points), and from the major axis (filled points) around 
the minor and major axes for the 71 galaxies ordered in right 
ascension considered in the case flmax = 3°. Bottom : Difference 
between number of QSOS around the minor and major axes. 
The continuum lines represent the +3cr level. 

is due to the real distribution, we have to carry out Monte Carlo 
simulations. 

4.1. Monte Carlo simulations with random position 
angle of the galaxies 

Given the present distribution of QSOs in Schneider et al. 
(2005), we have randomly generated the position angles of 
71 galaxies in the RC3 (the correlation of position angles, if 
any, are small and insignificant; see ^4.4L and we have mea- 
sured the same numbers Pk~s and £45. We have carried 50000 
Monte Carlo simulations of this type. The distribution of £45 
for the case without constraints is shown in Figure ^2 This 
distribution results in a Gaussian distribution with cr x 250. 
Our detection of anisotropy would be at the 3.5cr level. Table|3 
presents the results and the statistical significance for this and 
other cases with further constraints. 

It is remarkable that merely by removing the cases with 
z < 0.5 (only 9% of the cases) we increase the statistical signif- 
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Fig. 10. Top: Distribution of QSOs around NGC 3795. The line 
represents the minor axis direction. Bottom: Histogram repre- 
senting the number of QSOs as a function of position angle 
with respect to the minor axis of NGC 3795 within 3 deg of it. 
A constant dependence is excluded within the 7.2cr level. 



icance considerably (to 3.9cr). This is due to the negative signal 
of the anisotropy for these removed cases, as shown in Figure 
m As said, this better signal/noise of the anisotropy is also ob- 
served when comparing the two cases shown in Figure|8l 



Fig. 11. Results of the |£'45|, excess in the range 0^5° over 
the range 45-90° (with a total of 25 176 galaxy-QSO associa- 
tions), of a Monte Carlo simulation with random position an- 
gles for the RC3 galaxies. As shown, the real position angles 
give a significant anisotropy of £45 = 874, which is the 3.5cr 
level. 



4.2. The same test rejecting the three farthest outliers 

We repeat the statistical analysis rejecting the three farthest out- 
liers identified in Tableland Figure |5] By definition, outliers 
do not behave like the others, either because of the physics 
or because of certain uncontrolled factors; separating outliers 
might be a test to better understand the properties of the general 
population, and of the outliers themselves if they are physically 
different. We repeat the case only with the constraint Omax = 3°. 

We now remove three galaxies from Table |3 PGC 28741, 
PGC 36192 and PGC 51752, all of them with the highest pos- 
itive value of £45, so it is expected that we are going to reduce 
the signal of the anisotropy. The number of galaxies is now 68. 
The number of pairs is 23 970. £45 = 632, a = 0.097 ± 0.018, 
Pk-s = 4 X 10"^. The Monte Carlo simulations on the random 
orientations of the galaxies give a significance cr(PK-s) = 2.7, 
cr(£'45) - 2.6. It seems that there is nothing wrong with the 
previous statistics. The signal is still there in spite of the re- 
moval of the three farthest outliers. It is lower than 3.5cr and 
is reduced to 2.6 or 2.7cr, but this is precisely what would be 
expected if we remove members that are known to have the 
highest anisotropy. There is an average trend for a positive a 
that is not due to the presence of a single galaxy (or two or 
three) with some special circumstances. 
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Table 3. <t from the results of Monte Carlo simulations with random orientations of the galaxies for the measured values of Pk~s 
and |£'45| in the real distribution. 
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4.3. Tests with subsamples 

We now perform the same calculations with subsamples con- 
taining 90% of the identified QSO-galaxy pairs. If the sam- 
ple is pure, the anisotropy will decrease slightly owing to the 
smaller number of objects. If outliers have a strong effect, most 
of the simulations will have a slightly higher anisotropy and 
a few of them will have a much lower anisotropy (when the 
outliers are randomly rejected). 

For this test with 1000 different subsamples (each subsam- 
ple with 90% randomly selected pairs from the initial 25 176 
in 71 galaxies), we calculate the quantity E4s/Npai,-s, which 
gives a measure of the anisotropy. This time we do not cal- 
culate the significance for each subsample, because this is very 
time-consuming; in any case, the significance does not need 
to be calculated for each different subsample because it is 
more or less proportional to £45 /A^pairs ■ The results of these 
1000 subsamples are shown in Figure^l The median value of 
in thcsc 1000 subsamples is 0.0347, the same value 
that was obtained with the original total sample of 25 176 pairs. 
Nothing anomalous is present in this distribution of values of 
E45 /A^pahs SO we do not think that a few outliers are producing 
the anisotropy but that it is a characteristic of the average dis- 
tribution. All 1000 subsamples show a relative excess of QSOs 
with position angle less than 45 degrees of E45/N over 0.025. 

4.4. Monte Carlo simulations with random QSO 
distributions 

In the previous subsection, we showed that, given the actual 
distribution of SDSS QSOs, a random distribution of posi- 
tion angles in the RC3 galaxies should not show the level of 
anisotropy found; that is, the anisotropy is not due to the pe- 
culiarities of the QSO distribution itself. Now from the posi- 
tion and orientation of the 71 RC3 galaxies we generate a ran- 
dom Poissonian and homogeneous distribution of background 
sources through a Monte Carlo simulation and we shall see 
whether the observed level of anisotropy is attained. These 
should serve to test that the 7 1 RC3 galaxies have no signifi- 
cant bias responsible for the anisotropy. In fact, we have also 
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Fig. 12. Distribution of values of the relative excess of QSOs 
with position angle less than 45 degrees for 1000 subsam- 
ples (« from 1 to 1000), where each subsample has 90% ran- 
domly selected pairs from the initial 25 176 in 71 galaxies. The 
dashed line represents the value for the original 25 176 pairs: 

£45/A^pair,s = 0.0347. 



checked that the distribution of position angles in our sample 
(see Table O is compatible with being randomly distributed 
(according to the Kolmogorov-Smirnov test: a probability of 
83% for the hypothesis of null correlation). There is no signif- 
icant correlation of position angles, although, of course, some 
small correlation within the noise cannot be rejected; in any 
case, this small correlation would be negligible for our statisti- 
cal purposes. 

The results indicate that the probability of obtaining an 
anisotropy of the measured values is roughly equal to the prob- 
ability given by the Kolmogorov test. From 873 simulations, 
we got 89.9% with Pks > 0.1, 9.5% with 0.01 < Pk-s < 0.1 
and 0.5% with 0.001 < Pk-s < 0.01; none with Pk-s < 0.001. 
This distribution is to be expected if the probabilities obey the 
Kolmogorov-Smirnov distribution. These results indicate that 
the probability of getting the observed value Pk-s - 2x 10"^ 
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from a random distribution of QSOs is negligible. This means 
that there is no selection effect in the RC3 sample. 

4.5. Tests with other catalogues than Schneider et al. 
(2005) 

We have also carried out the same test with the spectroscopic 
galaxies of the 3rd SDSS release up to nig = 21 (333 314 
galaxies), instead of the QSO sample, and we have found that 
those galaxies with a difference of velocity with respect to RC3 
galaxies higher than 3000 km/s (to avoid their being associated 
with the main RC3 galaxy) have an insufficiently significant 
level of anisotropy. As always, the RC3 galaxies are spiral and 
nearly edge-on (inclination greater than 65 degrees). The re- 
sult gives a - -0.029 + 0.010 (more galaxies towards the ma- 
jor than towards the minor axis), Pk-s = 2 x 10 and the 
Monte Carlo simulations (with random orientation of the RC3- 
galaxies) show that this has a probability of 38% (equivalent to 
only 0.9O-). 

If instead of Schneider et al.'s (2005) vetted QSO catalogue 
we use the complete SDSS 3rd data release for QSOs (we also 
put the constraint that the difference in velocity of the QSO and 
galaxy must be higher than 3000 km/s), the anisotropy is de- 
tected with a = 0.1 1 1+0.022, Pk-s = 3 x 10"^, and, according 
to Monte Carlo simulations, the anisotropy has a level of 2.9cr. 
The reason for the lower significance of the anisotropy (with 
the Schneider et al. sample it was 3.5cr) is that, apart from fur- 
ther spurious contamination, it includes QSOs with M, > -22, 
which are mostly of low redshift, and this has negative a, as 
shown in Fig. |5j which reduces the average signal. With the 
added constraint z > 0.5, a = 0.154 + 0.024, Pk-s = 2 x 10"'' 
and the Monte Carlo simulations indicate anisotropy at a level 
of 3.6cr. 

We can also use the newest release, SDSS 4th (DR4, 
Adelman-McCarthy et al. 2006), which covers 4783 sq. de- 
grees and contains 61 049 QSOs with nig < 21.0. Although the 
area covered is only 14% greater than that covered by the DR3, 
the total number of RC3 galaxies that have circles of radius 3 
degrees totally covered is somewhat higher: 127 instead of 71, 
and the number of QSO-galaxy associations (some QSOs may 
have more than one association with a galaxy) is 47 271 instead 
25 176. This improves the statistics slightly. However, since the 
vetted catalogue has not yet been produced and will not be pro- 
duced until the 5th release is delivered (Schneider, priv. comm.) 
only the analysis of the anisotropy of the complete SDSS 4th 
data release is possible, as in the previous paragraph. The result 
is: a = 0.073±0.013,Pa:_s = 2 x 10"^ and according to Monte 
Carlo simulations the anisotropy has a level of 2.5cr. With the 
extra constraint z > 0.5, a = 0.106±0.012,P;f_s = 7x10"^ and 



the Monte Carlo simulations suggest anisotropy at a level 3.3cr. 
The average anisotropy is somewhat lower than with only DR3 
sources, possibly the first 71 galaxies have on average a higher 
anisotropy; in any case, the differences are roughly within the 
expected random fluctuations from one sample to the other 

Another search for anisotropy was made using the 2dF 
QSO spectroscopic catalogue (Croom et al. 2004) with nig < 
21 (23 660 QSOs in total). This catalogue has the disadvan- 
tage of a much lower coverage, around 700 square degrees, and 
its two strips have a width of 5 degrees, which does not allow 
complete circles of radius 3 degrees (we consider the full ar- 
eas of the two 75 deg x5 deg). We do the analysis with nearly 
edge-on spiral RC3 galaxies that have at least 80% of the cir- 
cles with 0max = 3° covered (26 galaxies) and we normalize 
the QSO counts by dividing by the covered area. This results 
in poorer statistics than SDSS, but it is illustrative that we also 
observe the anisotropy, although with very low significance. 
Figure [T3I shows a histogram of densities: a - 0.089 ± 0.022, 
£45 = 1.6 + 0.5 deg"^; a probability according to Monte Carlo 
simulations of 51% (0.7cr), which is too insignificant to be 
considered a detection of anisotropy but is compatible with 
the results from the SDSS survey. If instead of 80% minimum 
coverage, we use only galaxies with 90% minimum coverage, 
we have 13 galaxies, a = 0.109 ± 0.028, £45 = 2.1 ± 0.6 
deg"^; and a probability according to Monte Carlo simula- 
tions of 22% (1.2cr), still insignificant. For Smax = 2.5°, and 
a coverage greater than 90%: 20 galaxies, a - 0.084 ± 0.017, 
£45 = 1.2 + 0.6 deg"^; and a probability according to Monte 
Carlo simulations of 34% (l.Ocr), insignificant. 

5. Some considerations towards measuring the 
anisotropy 

5.1. Anisotropy as a function of linear distance 

One possible question about our measure is why we counted 
the QSOs within a fixed angular distance (3 degrees) instead 
of a fixed linear distance. Since we know the distance of the 
RC3 galaxies, we could count the QSOs within a linear pro- 
jected distance. Indeed, we tried to measure the anisotropy in 
such a way and we also detected the anisotropy but the statisti- 
cal significance was not as high as 3.5cr (3.9cr for zqso > 0.5). 
With identical criteria, using Schneider et al.'s (2005) QSO cat- 
alogues but for the maximum linear distance of 1 .8 Mpc (where 
we got the best values of the anisotropy) instead of 3 degrees 
of angular distance, the results are: 113 RC3-galaxies, 16572 
galaxy-QSO associations, a = 0.138 ±0.038, P/f-s = 2x10"^ 
with a significance according to Monte Carlo simulations of 
2.6cr. The two-dimensional distribution as a function of dis- 
tance is plotted in Figure fT4l For zqso > 0.5, there are 15 079 
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Fig. 13. Histogram representing the 2dF QSO density with 
niB < 21 within a circle of radius equal to 3 degrees of 26 RC3 
galaxies (more than 80% of the circle covered) as a function of 
position angle with respect to the minor axes of these galaxies. 



Fig. 14. Counts of QSOs (total: 16572) as a function of posi- 
tion with respect to the minor/major axis angles of the corre- 
sponding 113 galaxies for linear distances lower than 1.8 Mpc 
(see text). Counts were plotted in bins of 0. 12 Mpc xO. 12 Mpc 
(average around 1 10 QSOs per bin) and smoothed/interpolated. 
The clearest colours indicate higher density. 



galaxy/QSO associations, a = 0.172 + 0.042, Pks = 2 x 10"^ 
with a significance according to Monte Carlo simulations of 
3.1cr. 

This slightly lower significance of the anisotropy may be 
due to the following reasons: i) the selected number of very 
nearby galaxies (~ 10 Mpc) with the linear distance criterion 
is lower than the angular distance criterion because most of 
them do not follow the constraint of being totally covered by 
the SDSS 3rd release survey, a galaxy at ~ 10 Mpc would give 
a circle with radius ~ 9 degrees, which is very unlikely to be 
totally covered by the SDSS-3rd release survey (and, as said in 
^ there is a slight trend towards higher anisotropy for lower 
distances); ii) the number of QSOs per galaxy is lower for more 
distant galaxies, because they cover a lower area. These argu- 
ments can explain the numbers we obtained. 

One might restrict the range of distances to avoid this se- 
lection effect in some measure. For instance, if we take only the 
RC3 galaxies within the range of distance 20 < < 40 Mpc 
(around the median value of 32 Mpc; i.e. 0.005 < Zgai < 0.010 
instead of Zg-A < 0.050). We would then have 20 RC3 galax- 
ies, 9457 QSO-galaxy associations , Pks - 2x 10"^, Monte 
Carlo significance: 2.8-cr. For zqso > 0.5: 8631 QSO-galaxy 
associations, Pks = 1 x 10"^, Monte Carlo significance: 
3.3cr. The signal/noise is a little better with this restriction but 
still less than the values obtained with the condition < 3°. 
Here, the reason is possibly the lower number of galaxy-QSO 



associations used for the statistics. Whatever the reason, it 
seems that the condition 6 < 3° gets higher anisotropy than 
linear separation < 1.8 Mpc. 

5.2. Counting some QSOs twice or more times 

An interesting question is how to select QSO/galaxy pairs. At 
present, we have selected all the QSOs around a given RC3 
galaxy; consequently, if two galaxies are closer together than 
6 degrees, they will have some common QSOs within the cir- 
cle of radius 3 degrees; some QSOs count twice or more in 
the QSO/galaxy association. Since the position and orienta- 
tion of the galaxies are random, the different counts with the 
same QSO should be independent. Even in the case that they 
were not independent, the Monte Carlo simulations of ^.11 
would include the possible factors due to it, and the signifi- 
cance derived therein would contain them. We have estimated 
the anisotropics with other criteria in which a QSO counts 
only once: associating it with the closest/farthest galaxy in 
the case that two or more galaxies include it within their cir- 
cles, and the observed anisotropy has only slight differences: 
EAilNpair - 0.0328/0.0373 respectively for closcst/farthcst Cri- 
terion instead of 0.0347 (Table O. This is explicable perhaps 
because we would be adding a selection effect in which the 
QSOs are on average closer or farther from the galaxies; or by 
coordinates, in which a QSO is associated with the first/last 
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galaxy in the list of RC3, which is ordered by coordinates, and 
in this case we found more notable differences, which we at- 
tribute to the gradient of QSO density in right ascension present 
in SDSS survey. Summing up, the most objective criterion is 
to include all QSOs around a galaxy even if they are counted 
twice or more; other criteria can introduce extra bias; moreover, 
we have fewer counts and consequently less significance in the 
detection of the anisotropy, although whatever the criteria, the 
Monte Carlo simulations will include them and the significance 
will be correctly calculated. 

5.3. With incomplete circles 

So far, we have used circles that are covered completely by the 
SDSS survey (except for 2dF data). We have performed other 
measures with circles partially covered, normalizing the num- 
ber of counts with the area covered per segment, and we have 
observed that: i) it is a trivial result that the number of galaxies 
is higher for lower minimum allowed coverages; ii) lower min- 
imum coverages in the circles produce higher dispersion and 
lower average anisotropy in our case. When the non-covered 
holes are larger, random effects dominate [perhaps because we 
are taking border areas with significant gradients of complete- 
ness; also, the errors in the determination of the areas of the 
segments might contribute] and in our case the average random 
anisotropy is fortuitously negative. We have perfomed mea- 
sures only with galaxies of low coverage and we have tested 
that they produce a negative a, in which Monte Carlo simula- 
tions show that it is not significant. When we mix low-coverage 
with full-coverage circles we are wrecking/smoothing the ef- 
fect observed with the full-coverage circles alone. 

The best compromise between points i) and ii) to get the 
highest significance of the anisotropy was found for the con- 
straint of 98% as minimum allowed coverage. This gives 135 
galaxies (instead of 71 galaxies with total coverage), 48 131 
QSO-galaxy pairs , a = 0.133 ± 0.013, £45 = 0.837 deg^^ 
whose significance according to Monte Carlo simulations is 
4.4cr (a probability of 1 in 90 000). For z > 0.5 the significance 
according to Monte Carlo simulations is 4.8cr (a probability of 
1 in 600000). TableEJgives further values for other coverages. 

6. Interpretation 

This distribution of objects around galaxies is similar to that 
found by Zaritsky et al. (1997) and Azzaro et al. (2006) (the 
Holmberg effect), i.e. a higher density of objects in the mi- 
nor axis direction. However, their distribution was for satel- 
lite galaxies associated with the main galaxy and our objects 
are QSOs with very different redshifts with respect to the main 
galaxy. The projected scales are also different ~ 300-500 kpc 



for the satellites, and ~ 1-1.5 Mpc in our case. The Holmberg 
effect has been explained as a consequence of the preferential 
capture by the parent galaxy of satellites orbiting near the equa- 
torial plane. From a purely phenomenological point of view, 
the effect discovered here is similar to the Holmberg effect. 
Although we have not tried to build a detailed model, we have 
considered three possibilities: 

• Extinction: It might be that there is some extinction along 
the major axis of the galaxies that is high enough to reduce sig- 
nificantly the number of QSOs observed in that direction. The 
affected range of redshifts (z > 0.5) contain QSOs in which H^ 
lines cannot be detected in the optical SDSS survey (and Ly-or 
cannot be observed for z < 2.2), so lines like H„ (or Ly-a) are 
not used for their identification as QSOs when 0.5 < z < 2.2; 
instead, other, fainter lines are used, and they will be near the 
limit of detection at nig around 20. An excess of extinction 
along the major axis of ~ 0.1 mag could produce a selection 
effect in the SDSS survey that removes some of the objects to 
be classified as QSOs because the signal/noise of some lines 
becomes lower than the limit of detection. For example, if we 
have a QSO whose equivalent width of the most intense line 
has got a signal/noise of 4 without extinction, with extinction 
this signal/noise would be lower than 4, so the line of this ob- 
ject would not be catalogued and therefore this object would 
not be classified as a QSO. However, it is difficult to accept 
that some clouds along the major axis associated with a galaxy 
are as far as 1.5-2.0 Mpc (see Fig. I14> from it. Intergalactic 
clouds might reside at 1.5-2 Mpc from a galaxy but there is no 
reason in principle to think they are in the plane of the galac- 
tic disc rather than in a random distribution. A random distri- 
bution of intergalactic clouds would not produce anisotropy. 
Moreover, we have computed the polar and radial distribution 
of ^ - r colours of our sample of QSOs and we find no appre- 
ciable differences between the mean colour of objects at angles 
<> +45° from the minor axis of the galaxy (FigurelTslpresents 
the two-dimensional distribution of such colours; the average 
colour as a function of the position angle from the minor axis 
is {(g - r)) = 0.1727 + 0.0022 + 1.9 ± 4.3 x \Q-^PA{dsg)), so 
unless a grey-dust extinction is present, this explanation does 
not seem appropriate. 

• Gravitational lensing: As said in the introduction, previ- 
ous detection of a correlation between high redshift QSOs and 
low redshift galaxies has been tentatively attributed to effects 
of gravitational lensing in the halo of foreground galaxies. This 
solution does not work in general (Zhu et al. 1997; Tang & 
Zhang 2005), and it also has important problems in explaining 
the phenomenon of this paper: 

- The dependence of the anisotropy on the redshift of the 
QSOs (excess over z > 0.5 and defect over z < 0.5) can- 



14 



Lopez-Corredoira & Gutierrez: Anisotropy/QSOs 



Table 4. Statistics with different minimum allowed coverages of the circles with 6,^^^^ - 3°. 



Min. cover. 




N ■ 

^ ' pairs 


a 


£45 (deg 


o-{\E4s 


100% 


71 


25176 


0.132±0.017 


0.872 


3.5 


99% 


117 


41888 


0.124±0.015 


0.791 


4.0 


98% 


135 


48131 


0.133±0.013 


0.837 


4.4 


98% (zq > 0.5) 


135 


43839 


0.158±0.014 


0.903 


4.8 


97% 


143 


50938 


0.105±0.013 


0.677 


3.8 


95% 


162 


57000 


0.084±0.010 


0.546 


2.9 


90% 


188 


64791 


0.067±0.011 


0.452 


2.5 



Minor a 

3°- 




Fig. 15. Distribution of the average colour {(g-r)} as a function 
of position with respect to the minor/major axis angles of the 
corresponding 71 galaxies (see text). Average: 0.173. Bins of 
0.2° X 0.2°, smoothed/interpolated. The clearest colour indicate 
reddest colour. 

not be explained in terms of the distance of the source and 
the lens. All gravitational lens effects, whatever the model, 
have an Einstein radius 6e such that 6"^ is proportional to 
(e.g. n - 2 for microlensing, n — I for singular isother- 
mal sphere approximation), with D/s and Dj respectively 
the distance lens-source, source-observer. is nearly one 
for all QSOs at z > 0.08 (assuming a distance of the galaxy 
of 32 Mpc, ^ > 0.9). We cannot understand with this 
hypothesis the absence of anisotropy for most QSOs with 
z < 0.5. The dependence on the magnitude shown in Fig. 
|6lis also difficult to understand, we cannot understand why 
there is no lensing effect for the brightest QSOs. 
- The linear scales of the effect discovered here seem too 
large for the size of the halos as estimated in A^-body CDM 
or ACDM simulations. For instance, typical virial radii 
are ~ 100-300 kpc in these models. Although Prada et al. 



(2005) have demonstrated that actual halos extends at least 
up to 2-3 Ryir, the density at this distances is so low that 
it seems difficult to produce significant lensing effects. A 
point to note is that such simulations predict elongated ha- 
los (Allgood et al. 2006; Mandelbaum et al. 2006), whose 
major axes in the inner parts tend to be orthogonal to the 
galactic discs, while in the outer parts the halos tend to 
be randomly orientated. The anisotropy found in this work 
would need elongated structures perpendicular to the galac- 
tic discs and extending out to ~ 1.5 Mpc (see Fig. ll4> . 
- The amount of mass necessary is huge. Assuming as a first 
approximation that the increase in counts in SDSS is due 
to an increase in the magnitudes of the background sources 
(this is not exactly true, because the magnitudes of an object 
refer to the average flux of the object in a passband, while 
the SDSS detection criteria depend on the signal/noise of 
some lines, but the order of magnitude should not be very 
different), the average enhancement is {q) x 1 .07 (an in- 
crease of 7% in the counts) in the counts over ^ 14 deg^ 
(the region within 6 < 3° and polar angle within 45° from 
the minor axis). The magnification factor /i due to gravita- 
tional lensing [see, for instance, eq. (1) of Zhu et al. 1997] 
is related to q by means of 

A^(mA+2.51ogioyu) 1 
N{mt) n 

where is the cumulative QSO counts up to magnitude 
nth. With the Boyle et al. (2000) counts [A^ = 1981 - 
2\A.2mb + 5.792m^ (Lopez-Corredoira & Gutierrez 2004, 
appendix A); in the real distribution of QSOs, not in our 
sample which is incomplete], </i) ^ 1.28 for mt = 20.0 on 
average over ^ 14 deg^. For mt, - 19.5, it is (yu> =s 1.08 
and for mt - 20.5 there is no value of yu that gives an 
enhancement as high as ^ - 1.07. Let us consider the 
average value of mb ~ 20.0: (yu) ^ 1.28. With a single 
isothermal sphere approximation fi ^ "g^^ where 6e is 
the Einstein radius; whose average over a circle of radius 
^max >> S'e is ~ 1 + 1^. Therefore, in our case it 
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would be 0£ » 1100 arc seconds, which would require 
an excess of mass in the minor axis region with respect to 
the major axis region of a whole rich cluster of galaxies 
(for the typical galaxy distance of 32 Mpc) in a single spi- 
ral galaxy (Wu 1996). Even if the mass is distributed over 
many lenses, the total necessary mass is of the same order 
because {(ji - 1)> oc Y^iSej °^ Y^iMi - Mtotai, although the 
numbers can change depending on the size of the lenses. 
With a point-like microlensing approach (Paczynski 1986; 
Wu 1996), ^ ~ -^^==, M = O/Oe, whose average over 

a circle of radius Omax » Or is (/i) ~ I + it re- 
sults in 6e - 2900 arc seconds. In the case of microlens- 
ing, 0| is proportional to the mass (0| - for dis- 
tances of the QSOs much greater than the distance of the 
galaxy; Paczynski 1986; Wu 1996; Tang & Zhang 2005), 
<(yu - 1)) OC OC 2,.M,- = Mtotai, so again it does 
not matter whether the mass is distributed in one or many 
sources: the average magnification will be the same. In our 
case, we would need a mass M ~ 4 X lO'*" Mq for t/ = 32 
Mpc, an impossible value. 



• Non-cosmological redshifts: At least some QSOs are as- 
sociated with nearby parent galaxies. They would have abso- 
lute magnitudes between -10 and -14 in g. In this unortho- 
dox scenario, QSOs might be ejected with velocity enough to 
become gravitationally unbounded to the galaxy. Arp (1998a, 
Fig. 3-27) proposes ejection along the minor axis, although out 
to distances of 0.5 Mpc. If we were to assume this hypothesis 
to be correct, why do we find this effect at distances up to 2- 
3°? Maybe because our sample have relatively faint QSOs, and 
QSOs with small angular separation from a galaxy have higher 
luminosities on average (Dravskikh & Dravskikh 1996) that 
correspond to lower apparent magnitudes. Possibly the same 
effect of anisotropy would be observed for bright QSOs on 
scales less than 1°. Bright QSO excess might not be detected 
because the number of such QSOs in SDSS is very low for their 
excess to be detected statistically. Why do the QSOs with red- 
shift z > 0.5 preferably show this anisotropy and not other red- 
shifts? If Arp's hypothesis of ejection were right, maybe QSOs 
with lower redshifts are further away or turned into galaxies, 
but this places the maximum distance even beyond a: 1 .5 Mpc, 
so the problem with distance in Arp's hypothesis still remains. 
That interacting galaxies appear among the galaxies with the 
highest anisotropy in the distribution might point some relation 
between both phenomena. Therefore, Arp's hypothesis quali- 
tatively predicts the observed anisotropy but it fails in the dis- 
tance estimation of the QSOs from the galaxies by a factor of 
~ 3. 



Could the detected anisotropy be a statistical fluctuation? It 
might be, but it is a very low probability one: for instance 3.5cr 
means a probability 5 x 10 (with no preselection of QSOs 
and galaxies except the angle O^ilx - 3°). Moreover, Fig. |5] 
shows a clear dependence on redshift which is not expected for 
a random sample of QSOs (there is a clear non-random trend 
of anisotropy depending on z), so the significance increases to 
3.9cr only by removing the QSOs with z < 0.5 (9% of the to- 
tal). Is this possibly due to chance? We think that this is not 
very likely, but that is precisely the nub of the question. Here, 
we deliver this new challenge in the long-running puzzle of 
this old topic: the relationship between high redshift QSOs and 
nearby galaxies. In any case, since this is the first time that this 
anisotropy has been observed and the possibility of a chance 
fluctuation, although unlikely, is not totally discarded, we pre- 
fer at present to be prudent and just say that it is a "tentative" 
detection that should be corroborated by other groups before 
any extraordinary claims are made. 
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